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Abstract 

The renormalization of local dimension-4 operators containing a heavy and a light 
quark field at scales below the heavy-quark mass is discussed, using the formalism 
of the heavy-quark effective theory. The anomalous dimensions of these opera- 
tors and their mixing are calculated to two-loop order. Some phenomenological 
applications are briefly discussed. 



(Submitted to Physics Letters B) 



CERN-TH/97-304 
November 1997 



1 Introduction 



Local operators containing both heavy and hght quark fields exhibit an interesting be- 
haviour under renormalization at scales below the heavy-quark mass tuq. Large loga- 
rithms of the type asln(mQ//i) arise from the exchange of gluons that are "hard" with 
respect to the light quark but "soft" with respect to the heavy quark. Since such gluons 
see the heavy quark as a static colour source, the large logarithms can be summed to 
all orders in perturbation theory using an effective theory for static heavy quarks, the 
so-called heavy-quark effective theory (HQET) In the HQET, the 4-component 

heavy-quark field Q{x) is replaced by a velocity-dependent 2-component field hy{x) sat- 
isfying hy = hjj, where v is the velocity of the hadron containing the heavy quark. 
Because of the particular hierarchy of the mass scales involved, the renormalization of 
heavy-light operators in the HQET is called "hybrid" renormalization. Operators in the 
effective theory have a different evolution than in usual QCD. For instance, whereas the 
vector current q'y'^Q is conserved in QCD (i.e. its anomalous dimension vanishes), the 
corresponding current q'~f'^hy in the HQET has a nontrivial anomalous dimension @]-@], 
which governs the evolution for scales below the heavy-quark mass. It has been calcu- 
lated at the two-loop order in Refs. [0]-[pi. The renormalization of higher- dimensional 



heavy-light operators is known only at the leading logarithmic order [|TO|-|T^. Here we 
shall generalize these results and perform the calculation of the anomalous dimension 
matrix of local dimension-4 operators at the two-loop order. This calculation requires 
the evaluation of two-loop tensor integrals in the HQET that are infrared (IR) singular 
when one of the external lines is taken on-shell. A general algorithm to calculate these 
integrals has been developed recently [|l3|| . 

The matrix elements of local dimension-4 heavy-light operators appear at order l/mg 
in the heavy-quark expansion of heavy-meson decay constants [|l4l and weak transition 



form factors [|I^]. As such, they play an important role in many phenomenological 
applications of the HQET. The theoretical predictions for weak decay form factors involve 
operator matrix elements renormalized at the scale mq. Our results can then be used to 
rewrite these matrix elements in terms of ones renormalized at a scale /i -C mq, which 
may be identified with the scale at which a nonperturbative evaluation of these matrix 
elements is performed. Details of such applications will be discussed elsewhere. 



2 Operator renormalization 

We start by constructing, in the HQET, a basis of local dimension-4 operators containing 
a heavy and a light quark field. Since ultimately our interest is in the matrix elements of 
these operators between physical hadron states, it is sufficient to consider gauge- invariant 
operators that do not vanish by the equations of motions |T6|. Then, in the limit where 
the light-quark mass is set to zero, any dimension-4 operator must contain a covariant 
derivative acting on one of the quark fields. A basis of such operators, which closes under 
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renormalization, is provided byQ 



Oi = g r iDahy , 

O4 = q {iv- D^'ja'^T hy . 



(1) 



where F represents an arbitrary Dirac matrix, which may or may not contain the Lorentz 
index a. The Feynman rules of the HQET ensure that no Dirac matrices appear on the 
right-hand side of F. Moreover, for = only an even number of 7 matrices can 
appear on the left-hand side of F. 

The equations of motion, iv ■ D hy = and q {ipY = 0, imply that between physical 
states the operators (Oi — O2), O3 and O4 can be replaced by total derivatives of some 
lower-dimensional operators, which are renormalized multiplicatively and irrespective of 
their Dirac structure. The additional power of the external momentum carried by the 
operators does not affect the ultraviolet (UV) behaviour of loop diagrams. Moreover, the 
reparametrization invariance of the HQET |1^] ensures that also the multiplicative 
renormalization of the operator Oi is determined by the same Z factor that governs 
the renormalization of dimension-3 operators |1T[. However, there is a nontrivial mixing 
between Oi and the other three operators. 

We define a matrix Z of renormalization constants, which absorb the UV divergences 
in the matrix elements of the bare operators, by the relation Oi = J2j '^■ij bare- Then 
the matrix "y of the anomalous dimensions, which govern the scale dependence of the 
renormalized operators, is given by 



dZ 



7 



z-i 



dln/i 

From the arguments presented above, it follows that the matrix Z has the structure 



(2) 





Z2 


Z3 


Z4 







ZI + Z2 


Z3 


Z4 










Zi 







I 










) 



(3) 



where Zi is the renormalization constant of local dimension-3 currents, which is known 
at the two-loop order 0^0- Here we will calculate the remaining entries in the matrix 
Z with the same accuracy. 

In a minimal subtraction scheme, the renormalization constants are defined to remove 
the 1/e poles arising in the calculation of the bare Green functions with insertions of 
the operators Oi in dimensional regularization, i.e. in = 4 — 2e space-time dimensions. 



^We use the notation iD = id + ggA and {iD^ = —id + ggA. 
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Hence, 



A:=0 



The requirement that the anomalous dimension matrix be finite in the hmit e 



(4) 




imphes the relations [lIS 
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a. 



da. 



2a 



a. 



das 



da.. 



da. 



a. 



(5) 



where P{as) = dtts/dln/i^ is the f3 function. The first equation shows that the anomalous 
dimension matrix can be obtained from the coefficient of the 1/e pole in Z, whereas the 
second one implies a nontrivial constraint on the coefficient of the pole, arising at 
two-loop and higher order. 

Previous authors have calculated the renormalization constants Zi appearing in (0) 
at the one-loop order, finding that |]TD|, ITT 



Aire 



3 3 _^ _1 
"2'2'~ '~2 



Oial 



1,2,3,4. 



(6) 



This information is sufficient to reconstruct the matrix Z at order as- From the second 
relation in f|5|), it then follows that 



,(2) 



a. 



11 



6 



3 3 

'2 '2' 



-1, 



1 



0(a: 



(7) 



Here Ca = N, Cp = |(A^^ - l)/N and = | are the colour factors for an SU{N) 
gauge group, and Uf is the number of light-quark fiavours. This relation will provide a 
check on our two-loop results. 

The equations of motion impose two additional conditions on the renormalization 
constants. Substituting T = v°T' or F = 7"F' in the bare operator Oi in (||), and 
using that between physical states qV iv ■ Dhy = 0, whereas qipT'hv = ida [g7"F'/i^] is 
renormalized multiplicatively, we find that 



(Z-^)l2 + (Z-^)l3 + (Z 

:Z-i)i3-2(l-e)(Z-i)i4 



-1^ 



/14 



poles 



0, 



0, 



(8) 



where it is understood that in the second relation all pole terms proportional to 
must cancel. In terms of the constants Zj, we obtain the conditions 



0; A; > 1 



2Z 



(1) 



1 + 2Z['> + Z^2 +2^< 



(2) 



(9) 
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Figure 1: Two-loop diagrams contributing to the calculation of the 
renormalization constants Zi. The shaded circles represent one-loop 
insertions of the gluon self-energy. 



which will provide a highly nontrivial check on our two-loop results. Relations (§) allow 
to calculate two of the three constants Z2, and Z4 in terms of the third one. The fact 
that the calculation of one of these constants suffices to reconstruct the full matrix Z, 
using that Zi is fixed by reparametrization invariance, has been observed in Ref. 
however, the relations derived there are not fully correct beyond the one-loop order. 



3 Two-loop calculation 

To obtain the renormalization constants Zi at order a^, we calculate the insertions of the 
operator Oi into the amputated Green function with a heavy and a light quark to two- 
loop order. The relevant diagrams are shown in Figure We need to keep terms linear 
in the momentum p of the light quark, as only those contribute to the renormalization 
constants Z2, Z3 and Z4. However, for completeness we will also keep terms linear in 
the heavy-quark momentum ka, which contribute to Zi. In that way we will reproduce 
the known two-loop result for Zi as a check of our calculation. Because the pole parts 
are polynomial in the external momenta, we can first take a derivative with respect to 
p and k and then set p = 0, so that all integrals are of propagator type and depend 
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t + ka+ks+pt+p s + k 



Figure 2: Schematic representation of the R* operation. The black 
dots represent the original vertices, the dashed line the IR counterterm. 

on the single variable uo = v ■ k. However, this method fails for some of the diagrams, 
for which setting p = after differentiation leads to IR divergences. In these cases, we 
apply a variant of the so-called R* operation , which compensates these IR poles by 
a construction of counterterms for the IR divergent subgraphs. 

Consider, as an example, the first diagram in Figure |l[ Its contribution is proportional 
to the integral 

When linearizing this expression in p, we encounter an IR divergence from the region 
t — 0, which can be removed by adding and subtracting the IR subtraction term 

lIR _ fAd^,Ad+ (t + p)/3SaS^ 



The difference Di — D^^ can be evaluated using a naive linearization in p, since the 
behaviour of the IR subtraction term for t — is the same as that of the original 
integral. However, because the expansion of involves tadpole integrals that vanish 
in dimensional regularization, the difference {Di — -D{^)iinoarized coincides with the naive 
linearization of the original expression for Di. The contribution D^^, which is necessary 
to subtract the IR subdivergence of the original diagram, factorizes into an IR counter- 
term (the t integral) and the original diagram with the lines of the IR sensitive subgraph 
removed, and with t and p set to zero. The construction just described is schematically 
shown in Figure 0, where the dashed line represents the IR counterterm 

dH ^P{t + p)p_ ipj ^ ^,^,, r(rf/2)r(rf/2-i)r(2-^/2) ^ ^^2) 



{2t:Y f^i^t^pf (471)^/2^^^^ V{d-\ 
The remaining two-loop tensor integrals are of the general form 

\VS + Uj \vt + UjJ (-s2)"^(-t2)"4[_(3_^)2]-5 

= _7r'^(-2a;)2('^-"^-"4-"5)+"+™/;;l•••J;-(^;; {a,}) , 



/ 



( - s^f' ( - t^f [ - (s - pyf [ - (t - pyf [-{s- tyf' 
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Using the method of integration by parts 



, we obtain the recurrence relations 



(t/ _ «i - «3 - 2«5 + n) + «3 3+(4- - 5-) + ai 1+2-] {v; {a^}) 



(^d-ai-as- 2«5 + n) + «3 3+(4- - 5") + «i l+(2- - 5") Rl\:X{P\ 



(14) 



which allow us to express any two-loop integral in terms of degenerate integrals, which 
have ^2 = 0, 0^4 = or as = 0. Here l"*" is an operator raising the index ai by one 
unit etc., and [/ij] means that this index is omitted. The degenerate integrals can be 
related in a straightforward way to products of one-loop tensor integrals |13|. Using this 



technique, we have calculated the pole parts of the two-loop diagrams in the 't Hooft- 
Feynman gauge. The results are summarized in the first three columns of Table |l|. The 
renormalization scale /x is introduced by the replacement of the bare coupling constant 



-^bare 



with fL = /xeT-s/2(47r)" 



1/2 



with the renormalized one through the relation g 
in the MS scheme. 

The two-loop diagrams in Figure |l] contain subdivergences, which must be subtracted 
by UV counterterms. In addition to the one-loop counterterms for the quark and gluon 
propagators and vertices, local operator counterterms are required. To find these, we 
calculate at the one-loop order all insertions of the operator Oi into the amputated Green 
functions with a non-negative degree of divergence. In our case, those are the two- and 
three-point functions with field content qh and qhA. We find that in the 't Hooft- 
Feynman gauge the UV divergences of these functions are removed by the counterterms 



a 



Oi - -O2 + O3 + -O^ -Ca — 
47re V 2 ^ 2 7 47re 



(15) 



Since the two-loop calculation was performed off-shell, a gauge-dependent operator has 
to be included in addition to the operators Oi introduced in (|l|). However, there is no 
need to include operators that vanish by the equations of motion [|T^]. The results for the 
sum of all counterterm contributions are summarized in the last two columns of Table |I[ 
When these contributions are added to the result for the sum of the two-loop diagrams, 
all nonlocal 1/e divergences proportional to ln(— 2co'//i) and ln(— p^//i^) cancel. This is 
a nontrivial check of our calculation. 

The sum of the two-loop diagrams plus their counterterms determines the two-loop 
coefficients in the products —Z]^'^Zy^Zi. To obtain the results for the renormalization 
constants at order a^, we have to account for wave-function renormalization of the ex- 
ternal quark fields. Using the one-loop expressions for the Zi factors in (P), as well as 
the known one- and two-loop (for Zi) wave-function renormalization constants of the 
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Table 1: Sum of two-loop and counterterm contributions in units of (ag/ Air)'^ 



Structure 












Pa 






1 1 

2e2 e 


J_ + 9. 

2e2 ^ 4e 


5 
4e2 


1 1 
2e2 e 






3e2 V 36 ^ 12 y e 


1 1 

2e2 € 


19 151 
12e2 18e 


11 

6e2 


11 I 22 
3e2 + 3e 




( n / n T) c 
P ± J? 1 t f 


3e2 ^ 18e 





2 _|_ 31 

3e2 9e 


2 

3e2 


4 8 
3e2 3e 


V-pVa 




5 1 / 17 1 27r2 \ 1 
4e2 Vl2 9 J e 


1 5 





3_ , J_ 

2e2 2e 


£2 ^ e 






il + f_Z _ zl!") 1 

6e2 ^ V 9 18 / e 








11 1 11 
3e2 3e 







C pTp n f 


2,5 
3e2 9e 








4 4 
3e2 3e 







CI 


8e2 ^ 9 J e 


1 3 

2e2 e 





3_ , J_ 

4e2 ^ 2e 


2e2 ^ e 




CpCa 


11 _|_ / 19 7r2\ 1 
12e2 """1,36 36 / e 








11 

6e2 







CpTpUf 


1 1 

3e2 18e 








2 

3e2 







CI 










1 

e2 







CpCa 










2 

e2 







CpTpUf 


















quark fields |22, 23], we obtain the final results: 



SCpGs 

Sire 



An 



5 An' 



^ / 11 49 7r2\ ^ /2 5 

+ C'a \ ] +Tpnf{ 

I 6e 36 9 / n3e 9 



8ne 



An 



' f 3 5 27r2^ 



47re 



1 + 



47r 

11 26 vr^ 



3 13 27r2^ 
4^ ~ 12 ^ 
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6e 9 



, +Tpnf , 
18/ n3e 9 



87re 



1 + 



47r 



4e 



1 
3 



27r2 



^ / 11 41 TT^ ^ /2 5 



7 



^ / 11 19 7r2\ ^ /2 1 



(16) 



Our result for agrees with that derived in Refs. The two-loop results for 

the remaining constants Zi are new. Note that these expressions do indeed obey the 
relations in (0) and (^. This a highly nontrivial check, which gives us confidence in the 
correctness of our results. 

The anomalous dimension matrix 7 governing the scale dependence and mixing of 
the renormalized operators can be obtained using the first relation in (^). Since the 
matrix coefficients in the perturbative expansion 



have the same texture as the matrix Z in @, it is sufficient to quote the entries of the 
first row, which we denote as 7^ with i = 1,2, 3, 4. We obtain 

(70)1, = 7,0 = (-3, 3, -2,-1) , 



57r 



2 


47r2\ 


'3 


9 ) 


104 




9 




+ 9 



^2/5 Svr^ 47r2 13 

(71)1, = 7i,i = '^^l2~~'~ ~'y 

/ 49 2n^ 41 104 2n^ 19 vr^' 

/lO 10 28 2\ 

We note that 72 + 73 + 74 = as a consequence of the first relation in (||) , and that the 
one-loop matrix coefficient 70 satisfies the simple relation 70^ = — 3Cf7o. 



4 Nontrivial basis transformations 

The results derived in this work are sufficient to calculate, at the two-loop order, the 
hybrid renormalization of any local dimension-4 operator containing a heavy and a light 
quark field. However, it some cases the choice of the operator basis in (|lD may not be 
the most convenient one. Whereas the transformation between one operator basis and 
another is trivial at the one-loop order, it may become subtle at the next-to-leading order, 
because in dimensional regularization the relations between the operators of different 
bases may depend on e. 

Consider the general case where the operator basis {Oj} is replaced by a new basis 
{Qj} with a linear relation of the form Oi = J2j ^ij{^)Qj- Depending on the choice 
of the Dirac matrix F, some of the operators Oi may not be independent, and thus it 
may happen that the new basis contains less than four operators. Therefore, in general 
the transformation matrix R(e) is a 4 x n matrix with n < 4. We introduce an n x 4 
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left-inverse of this matrix such that L(e) R(e) = 1. It then follows that the nxn matrix 
2 of renormalization constants in the new operator basis satisfies the relation 



r-l 



1 + L(e) Z 



r-l 



1 R(e) 



poles 



(19) 



where only pole terms are kept on the right-hand side. The anomalous dimension ma- 
trix in the new basis is determined by the terms of order 1/e in Z. Expanding the 
transformation matrices as R(e) = J2n Rn^" and L(e) = J2n Ln^", "we obtain 



^(1) = Lo Z(i) Ro + Li[z(2) - (Z(i))2l Ro + Lofz^^) - {Z^^)f] + . . . 



(20) 



where the ellipses represent terms that do not contribute at the two-loop order. Denoting 
by 7 the anomalous dimension matrix in the new operator basis, we obtain, using (^, 



70 = Lo 7o Ro , 

71 = Lo 7i Ro - Li(|7o^ + PoJo) Ro - Lo(|7o^ + Polo) Ri • 



(21) 



Alternatively, we may evaluate ([T9|) directly to get the inverse matrix and then use 
the first relation in (|^) to calculate the anomalous dimension matrix. 

As an example, we discuss the mixing of the local dimension-4 operators appearing 
at order l/mg in the heavy-quark expansion of the vector current| = q^^Q- There 



are six such operators, which can be chosen in the form [|T0|, pT 

Qi = qY^phv, QA = q{iv-DyYh 

Q2 = qvnpK, Q^ = q{ivD)^v^h 

Q^ = qiD^K, Q^ = q{tD^yK. 



(22) 



Each of the first three operators has the structure of Oi, with the substitutions T = 7^7°^, 
f'^7° and g'^", respectively. In all three cases, the corresponding operators O2, O3 and 
O4 can be written as hnear combinations of Q4, Q5 and Qq. For instance, the basis 
transformation for the case F = 7'^7" is 



(23) 



Using (|T^) and the constraints (|]) imposed by the equations of motion, we find for the 
corresponding matrix of renormalization constants 







/ 1 \ 




(QA 


02 




2 




Qa 






1 










\ 2e 4(1 -e) 







4,5,6) 



/ 2Z4 2Z3 2Z2 \ 

Zi 

Zi 

V Z4 Z3 Z1 + Z2J 



(24) 



^The discussion for the axial vector current is identical, provided the basis operators are chosen in a 
convenient way 
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where the subscript indicates the operators Qj whose renormahzation is accomphshed 
by this matrix. Similarly, for the other two cases F = and g'^"' we find, respectively, 



'(2,5) 




'(3,4,5,6) 






V 



Zi 


^4 



Zs 


Zi 

Z3 



Z2 






(25) 



Note that the 3x3 submatrix renormalizing the operators Q4, and Qq is the same 
in ^(1,4,5,6) and .2(34 5 6). This provides a nontrivial check of the formalism. The above 
results can thus be combined into a 6 x 6 matrix of renormahzation constants for the six 



vector- current operators in (^2]). The corresponding anomalous dimension matrix 7v 
has the block form 

A 



7v = 7i 1 



where 



/ 274 273 272 


V 74 73 72 



B 



B 




(26) 



(27) 



The exact two-loop expressions for these matrices, which are obtained from (|18|), are an 
important result of our work. 



5 Conclusion 

We have calculated, at the two-loop order in the HQET, the renormahzation and mixing 
of local dimension-4 operators containing a heavy and a light quark field. Our formalism 
allows for an arbitrary Dirac and Lorentz structure of these operators. The relevant two- 
loop diagrams have been evaluated with the help of tensor recurrence relations, once the 
IR divergent subgraphs have been regulated using IR counterterms. We find that a set 
of (up to) four operators closes under renormahzation. The equations of motion imply 
nontrivial relations between the anomalous dimensions of these operators, which we have 
derived. We have also shown how general basis transformations, which may depend on 
the dimensional regulator e, can be implemented in our approach. 

The matrix elements of local dimension-4 heavy-light operators play an important 
role in many phenomenological applications of the heavy-quark expansion. In partic- 
ular, they appear at order l/mg in the heavy-quark expansion of heavy- meson decay 
constants, and of the semileptonic form factors describing, e.g., 13 ^ it £u and B ^ piu 
transitions. Our results are an important step towards a full next-to-leading order anal- 
ysis of these quantities. However, still missing is the two-loop mixing of some nonlocal 



operators with the local operators considered here. This will be discussed elsewhere [24 
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